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INTEGRAL  EQUATION  OP  THE  THEORY 
OF  LIFTING  SURFACES 

N.  N.  Polyak  ho v 


This  work,  which  is  a continuation  and  a generalization  of 
work  [1],  sets  up  a new  form  of  an  integral  equation  of  the 
lifting  surface  and  examines  the  method  of  its  approximate  solution. 

§ 1.  Setting  up  the  basic  equation.  Let  us  assume  that  we 
have  a supporting  surface  which  represents  a set  of  rectangular 
lines  lying  in  the  plane  y=0  and  which  are  covered  by  a continuous- 
ly distributed  system  of  the  /7-shaped  vortices  with  density  y.  We  as- 
sume that  the  length  of  each  line  along  the  x axis  is  equal  to  2 cl, 
the  width  is  equal  to  dz , and  the  span  along  the  z axis  is  2L. 


The  inductive  velocity  caused  by  one  /T-shaped  vortex,  located 
as  shown  in  the  figure,  at  point  N'(x',  0,  z')  which  lies  in  the 
middle  of  the  line  with  number  k is  equal  to  [1] 
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i+nr*). 


Let  us  replace  the  variable  assuming  that 


+ a (*)+&*(#, 

,.*• -T4-  v * • • . v 


then 


p(*H^(*><  i +t)+v-f4 


where  x^  is  the  coordinate  of  the  leading  edge  in  the  z segment, 
while  the  dashes  denote  the  values  taken  in  the  £'  segment.  It  is 
evident  that  p characterizes  the  relative  shift,  while  q character- 
izes the  relative  narrowing  of  the  lines.  The  vortices  of  all  lines 
with  the  exeption  of  the  line  with  number  k inside  of  which  lies 
the  point  N' , will  yield  the  velocity 


where  all  z pertain  to  the  L value  and 


Vortices  of  the  line  with  number  k will  elicit  the  velocity 


Let  us  use  the  approximation  equality 


[ y \ — z1  dz 

J i*-*y 


+ 


which  will  be  more  accurate  the  smaller  the  «.  Then  we  will  obtain 


As  a result 
obtain 


— *>-^4 


^~r,  + 'f'  . <!+/')+ 
r=-u®  /l-- 


1 

-1  s*+ 1 


Vi  -*« 

isolating  the  peculiarity  and  integrating,  we  will 


(4) 


where 


is  a regular  function  since  1 — Kt  disappears  on  the  order  of 
(S'  — I)2  when  6-»S'. 

The  integrals  which  are  in  the  first  brackets  can  be  determined 
numerically  with  given  t , z'  . The  integrals  which  are  in  the 
second  brackets  can  be  determined  only  after  the  values  of  y in 
all  sections  can  be  determined  in  the  first  approximation. 

On  the  basis  of  (4)  the  equation  of  impermeability  at  the  point 
N'  will  assume  the  form  of  the  following  integral  equation: 


+i 


Vi 


(5) 


where  F is  a given  function.  The  equation  of  this  form  was  obtained 
earlier,  but  now  its  kernel  K has  a simpler  form,  which  is  important 
for  obtaining  the  approximate  solution  (see  [1]).  The  integral 
equation  (5)  permits  an  innumerable  number  of  solutions*.  From  this 
number  it  is  important  to  isolate  a solution  which  corresponds  to 
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the  smooth  flow  near  the  trailing  edge,  which  leads  to  the  require- 
ment that  the  density  Y*  disappears  on  this  edge.  The  disappearance 
should  have  a strictly  defined  order,  since  the  principle  value  of 
the  integral  which  enters  the  left  side  of  equation  (5)  exists  only 
in  certain  classes  of  functions,  in  which  the  integral  permits  the 
inversion.  The  inversion  formulas  have  the  following  form  (see  [2] 

§ 88) 

a)  ,.<0=^!4  SY&-&: 

-1 

b)  ; 

-t 

where  the  right  side  of  (5)  is  expressed  in  terms  of  <l>. 

After  the  inversion  of  the  integral  of  equation  (5)  under  the 
condition  that  the  postulate  of  S.  A.  Chaplygin  (*i(4l)*0)  is  satis- 
fied according  to  formula  a),  we  obtain 

T .(0-1.-  — i)— s-Vr=v  \ «<*>♦«.  '•  d*.  (6) 

2*1'  V l — r ' 1\ 

where  V.  corresponds  to  an  ifinite  elongation  and 

+i  

f l/l  ^ V K(i,  V.  z')<K' 


. 1 i/i  - i ' K (t,  V.  z')  & 

■ ~ ,i  V i - i'  — / 


One  can  attempt  to  solve  this  equation  by  the  iterations  method, 
having  assumed  that  in  the  expression  for  K the  Y”*V*  in  the  first 
approximation.  Then  equation  (6)  will  transform  into  the  Fredholm 
equation  of  the  second  order  with  the  regular  kernel.  In  order  to 
solve  it  by  the  iterations  method,  it  is  convenient  to  present  the 
solution  in  the  following  form 

T.(*)-T»<«>+  Vt=1+  *1*  , 

where  A-yJ,(±!)  = 0. 

After  the  substitution,  the  equation  (6)  will  assume  the  form 
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(7) 


7.(0 


T,»  ~ ^T«> 


: + HPy^rJM 


where 


+i 


At;,  o-  ?'  VrH** 

-t 


When  solving  equation  (7)  by  the  iterations  method  one  should  as- 
sume that  is  equal  to  (AT.*=0)  in  the  zero  aonroximation. 

This  enables  one  to  determine  the  in  all  sections  z'k  in  the 

first  approximation,  i.  e.  , to  determine  *,(S,  z),  then  find  the 
differences  l(t,  *)  — 7*($t  x)  and  Art*  and  proceed  with  the  calcula- 

tion of  the  second  approximation.  The  solution  becomes  completely 
elementary  if  one  uses  the  Vayzinger  [Weisinger]  hyDothesis  and  one 
assumes  that  £ is  equal  to  unity  when  calculating  J-|  and  A?  and 
assume  that  I'  is  equal  — 0.5.  From  formula  (7)  it  is  evident  that 
T (A  t')  vanishes  just  like  (1  - z,2)°*^  when  z'->±l,  and  the 
derivative  goes  to  infinity  just  as  (1  - z'2)-0’^  when  z'-*  ±1. 

We  can  prove  to  ourselves  that  these  boundary  conditions  are  a 
direct  consequence  of  the  vortex  method.  Actually  the  vortex  sheet, 
which  corresponds  to  the  curve  £=const,  by  intersecting  with  the 
Treffts  [sp.  unconfirmed]  plane  at  infinity  which  is  perpendicular 
to  it  yields  a rectilinear  segment  with  the  span  (-1,  1)  which  is 
covered  by  a vortex  layer  with  density  . The  inductive  velocity 
at  the  points  of  this  segment  is  equal  to 


djL  4a 
da  s’  — * 


Since  due  to  symmetry  the  circulation  near  the  examined  segment  is 
zero,  the  interval  should  be  inverted  according  to  formula  c)  which 
in  this  case  will  have  the  form 
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which  corresponds  to  the  boundary  conditions  formulated  above.  From 
the  aforesaid  it  is  evident  that  the  boundary  conditions  for  y , both 
at  the  leading  edge  and  the  trailing  edge  and  also  on  lateral  edges, 
are  the  consequence  of  the  introduction  of  the  vortex  system  for 
which  the  appearance  of  singular  intervals  which  permit  strictly 
defined  conditions  at  the  boundaries  of  the  integration  region  is 
characteristic. 

§ 2 . Substituting  the  Integral  equation  with  a system  of  al- 
gebraic equations.  This  operation  should  be  carried  out,  understand- 
ably, while  adhering  to  the  boundary  conditions  indicated  above. 

The  presence  of  the  singular  integral  in  equation  (5)  forces  one 
to  seek  the  solution  in  the  form 

U=Ao  Krn  + I*’  (8) 


where  *(±n  should  vanish  just  as  vr=r  when  5 — ► ± 1 . Sub- 

stituting (8)  in  (5)  and  eliminating  the  subscript  k,  we  obtain 


(9) 


where  K «K  (1,5'.  z') : 2k'  ( 1 - 

To  calculate  the  last  integrals  in  the  left  and  right  sides 
of  equation  (9)  we  divide  the  interval  of  integration  into  N partial 
intervals  (— 1,  5t),  (5t,  5,) . . . (5,v_i » ')  . We  will  assume  that  the  function 
K (5,5',  z’ ) is  known  in  each  partial  interval,  having  assumed  that 
in  the  first  approximation  -{*T*  • In  the  first  and  last  intervals 
we  approximate  n according  to  the  parabolic  law: 

% = • V — V-I  y ’ 


and  in  the  remaining  - according  to  the  linear  law: 


’L-’W-.+Ci-V.jljUfe. 


where  y - values  of  M at  the  points  of  division.  In  this  approxima- 
tion the  second  interval  on  the  left  side  of  equation  (9)  assumes 
the  form 


2b 


+ 1 ' Af-l  ' JH 


N- 1 


where  C*m  represented  by  the  integrals  are  calculated  very  simply5 
analytically  for  any  including  the  interval  containing  the 

singularity.  The  integral  on  the  right  side  is  brought  to  the  fol- 
lowing form  exactly  in  the  same  way 

+>  *-i 

*«»-i  *-• 


The  integrals  which  enter  the  C"m  can  be  calculated  using  the  for- 
mulas of  mechanical  quadratures,  since  they  do  not  have  peculiari- 
ties. As  a result,  equation  (9)  will  proceeed  to  the  following 
system  of  algebraic  equations: 

1 N- 1 

a,c0(?\  *')+  2 \c^v:  -^or.  n ( io ) 

In  this  case  the  ^o.  "fo,  % are  unknown.  Their  number  is 

equal  to  the  number  of  intervals  and,  therefore,  if  inside  each  of 
these  we  take  a point  with  the  coordinate  we  will  have  a system 
which  will  enable  us  to  determine  all  the  unknowns.  The  density  y 
is  determined  by  formula  (8).  This  approximate  solution  will  satis- 
fy the  required  boundary  conditions  of  the  vortex  method. 

§ 3*  Method  of  discrete  vortices.  There  are  two  concepts  regard- 
ing this  method.  One  of  these  is  the  Folkner  concept  and  it  entails 
the  fact  that  the  supporting  surface  is  broken  down  into  a series 
of  lines  which  are  parallel  to  the  direction  of  motion.  Each  of  the 
lines  is  covered  by  a finite  number  of  H- like  vortices  with  the 

circulations  f~ , where  m - line  number,  n - point  number  in  this 

mn  ’ 
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line.  In  each  section  the  values  of/"  join  the  coefficients  of 

' mn  " 

expansion  of  the  density  v*  in  the  series 

^A. slnd>’  (n) 

where 

5 — '—Co®  I,  — A, 9 ■j"  Ait^n  4*  4"  > • • * * — 0,  1.2,  * • * • 

The  connection  between  the  and  is  established  by  equating 

the  total  circulations  originating  from  the  discrete  and  continuous 
distributions  in  each  section.  The  same  thing  is  done  with  the 
inductive  velocities  at  the  assigned  points.  In  this  way  all  /"mn 
can  be  expressed  in  terms  of  A„*  and,  consequently,  total  inductive 
velocity.  Folkner  does  not  write  the  integral  eauation  explicitly, 
and  for  determining  the  A**  he  uses  the  system  of  equations  which 
expresses  the  condition  of  impermeability  at  the  various  points  of 
the  sections  used.  The  boundary  conditions  which  correspond  to  the 
vortex  method  are  satisfied  in  this  method. 

Another  approach  towards  the  method  of  discrete  vortices  con- 
sists of  the  fact  that  the  expansion  (II)  is  not  introduced  and  the 
circulations  of  /l-shaped  vortices  (the  usual  rectangular  or  skewed) 
remains  under  the  conditions  of  impermeability  so  that  what  we  get 
is  a system  of  algebraic  equations  not  relative  to  A,M,  as  Folkner 
has,  but  relative  to  r.  . This  system  has  the  form 

m 

(12) 

<-i 

where  i - vortex  number,  j - number  of  a point,  in  this  case  the 
numeration  proceeds  successively.  Furthermore,  segments  of  lifting 
vortices  are  located  at  0.25h,  and  the  condition  of  impermeability 
is  fulfilled  at  points  0.75h,  where  h - length  of  partial  inter- 
vals into  which  the  cross-sections  z=const  are  broken  down. 

Using  the  example  of  a two-dimensional  stationary  flow  which 
will  be  described  by  equation  (5)  with  the  zero  value  of  the  integral 
on  the  right  side  of  (A=*ooj,  we  can  see  that  system  (12)  cannot 
give  us  a correct  solution  for  the  problem.  Actually,  this  system 
obtained  from  an  integral  equation  will  have  the  form 
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(13) 


* 'i  « 

2 *‘J  f V*  = Fj'  Kij  = {t'i-0~  • 

Firstly,  the  substitution  of  the  kernel  using  formula  (13)  in  the 
interval  where  it  becomes  unlimited  (£J  - coordinates  of  discrete 

vortices)  is  totally  inadmissible  and,  secondly,  even  we  accept 
structure  (12)  for  the  initial  system  the  theory  of  integral  equa- 
tions says  the  following  concerning  the  correspondance  between  (12) 
and  (13)  (see  [3],  pg.  12):  " even  if  the  determinant  of  system 

(12)  does  not  equal  zero  and  the  system  has  one  solution,  system 

(13)  will  have  an  infinitely  large  number  of  solutions,  since  only 

the  average  values  of  the  y function  in  the  intervals  (—1,  |i)  , 

(6..  6»)  , etc.  are  determined  uniquely". 

It  is  known  that  in  order  to  select  a particular  solution  from 
an  infinite  number  of  solutions  which  would  correspond  to  the 
physical  conditions,  it  is  necessary  to  put  forth  the  boundary 
conditions  which  would  correspond  to  the  vortex  method.  Their 
satisfaction  must  be  guaranteed,  but  system  (12)  which  does  not 
reflect  to  any  degree  the  conditions  at  the  boundaries  does  not 
fulfil  this.  It  is  necessary  to  remember  that  the  approximating 
function  must  satisfy  the  same  boundary  conditions  as  the  approxima- 
tion function.  It  is  assumed  that  the  vanishing  of  the  vortex 
density  at  the  trailing  edge  in  the  limit  will  be  fulfilled  auto- 
matically if  the  last  oontrol  points  in  each  section  are  taken 
between  the  trailing  edge  and  the  last  discrete  vortex.  There  are 
no  substantiations  for  this  assertion.  In  this  case  it  is  necessary 
to  remember  that  what  is  important  here  is  not  just  the  vanishing, 
but  the  vanishing  of  the  rigidly  defined  order.  The  selection  of 
the  last  point  does  not  correspond  to  the  structure  of  a particular 
integral  and,  therefore,  cannot  supply  the  necessary  order  of 
disappearance.  Actually,  this  integral  exists  only  in  the  concept 
of  the  principal  value  which  is  obtained  for  the  point  lying  bet- 
ween the  two  segments  of  the  vortex  layer  to  which  corresponds  the 
point  lying  between  two  vortices  in  the  discrete  scheme.  The  "last" 
point  does  not  satisfy  this  condition.  Thus,  system  (12)  cannot 
replace  the  initial  integral  equation  with  the  boundary  conditions 
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which  correspond  to  it.  Therefore,  the  assertion  that  the  method 
of  discrete  vortices  in  the  theory  of  the  lifting  surface  " leads, 
in  all  cases,  to  the  substitution  of  two-  or  three-dimentional 
integrodif ferentia]  equations  with  the  systems  of  algebraic  equa- 
tions" (see  [4],  pg.  222)  is  not  true.  Such  a substitution  for  the 
equations  of  the  first  order  without  the  guaranteed  satisfaction  of 
the  rigidly  defined  boundary  conditions  is  not  lawful.  The  solution 
depends  on  the  nature  of  distribution  of  the  vortices  at  the  cal- 
culated points  which  is  assigned  quite  arbitrarily. 
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Summary  S 

The  paper  represents  the  simplification  of  the  article  [»].  The  integral  equation  accor- 
ding to  lilting  surface  theory  has  a form  (5)  which  can  be  transformed  to  the  form  (7). 
It  is  possible  to  solve  (71  by  iteration.  The  second  paragraph  gives  a system  of  linear 
equations  (10)  which  is  approximately  equivalent  to  integral  equation  (5).  The  vorticity 
distribution  per  unit  chord  y is  given  by  (8).  q(±l.  e)-0.  ' 
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